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I. INTRODUCTION 

Understanding the quark-gluon plasma created in nucleus-nucleus collisions 0, 0, 01 at the Relativistic Heavy Ion 
Collider (RHIC) is a challenging task. Many measurements have been carried out to provide insight on the properties 
of that dense QCD matter, yet a few years after the discovery, the question whether the plasma formed at RHIC is 
weakly or strongly coupled remains unanswered. On the one hand, results on bulk quantities like the elliptic flow 
vi or the shear viscosity to entropy ratio 0, 0j @| fit within a perfect- liquid picture. The large flow and 
low viscosity are naturally understood from strong cross-sections. This quickly prompted claims that the plasma is 
strongly coupled, although from a broader perspective, these issues are still being debated. On the other hand, one 
would like to establish a global picture and be able to also answer the question using hard probes. 

Hard probes are believed to be ideal processes to study the properties of the quark-gluon plasma, they are thought 
to be understood well enough to provide clean measurements. Observables built to measure medium effects on particle 
production, like the nuclear modification factors Raa, are not easily reproduced in perturbative QCD (pQCD), even 
for hard production. For instance to reproduce light-hadron Raa^ within the pQCD framework of medium- induced 
energy loss [1, [l^, [HI [13, the value of the so-called jet quenching parameter q is adjusted to 5 — 10 GeV^/fm 
[1^ . It is otherwise estimated to be smaller (1 — 3 GeV^/fm) for a weakly coupled plasma at RHIC temperatures. 
Results on heavy-quark energy loss are also puzzling, high— p_L electrons from charm and bottom mesons decays seem 
to indicate a similar suppression for light, charm and bottom quarks. By contrast in pQCD, the heavier the quark the 
weaker the suppression. While our current understanding of the pQCD picture must be still be improved (for instance 
collisional energy loss should be included as well), it is unclear if a pQCD approach can describe the suppression of 
high— particles at RHIC. Recent review and extensive discussion can be found in Refs. [Tsl. [l6l. \v\ . 

This motivates to think about strongly coupled plasmas. Addressing the strong coupling dynamics in QCD is 
an outstanding problem, the available tools are quite limited. However for the A/" = 4 Super Yang Mills (SYM) 
theory, the AdS/CFT correspondence [l^, [l^, [l^l is a powerful approach. Essentially, in the large— A^c hmit, the 
correspondence maps the quantum dynamics of the gauge theory at strong coupling into classical gravity dynamics 
in the fifth dimension of a curved space-time. While the SYM theory is quite different from QCD (it is highly 
supersymetric and conformal), let us assume that the SYM plasma approximates well the QCD plasma, just above 
Tc, the temperature of the phase transiton. In fact it has been argued that some SYM results hold for any gauge 
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theory with a gravity dual, hke the lower bound of the shear viscosity [1, 0|- If the plasma formed at RHIC really 
is strongly-coupled, this insight gained in AdS/CFT calculations like the one in this paper could help developing a 
qualitative understanding. 

Recently, there have been a lot of AdS/CFT studies on heavy quark propagating in the SY M p lasma. The energy 
loss is computed in Refs. [2l|, [2p. the momentum broadening is calculated in Refs. [H, 113, l25j . the j e t q uenching 
parameter is evaluated in Ref. [26l [27l . [28l | and the stress-tensor of the system is derived in Ref . [291 IsoL [biI. Is^ . 
Moreover, the saturation momentum and structure functions of SYM plasma are studied in Ref. [33l. l34l ISa. l36j and 
jet fragmentation and energy correlation are investigated in Ref. 37, 38, 39, 40]. Finally, there have been a few 
interesting studies which com par e the predictions of heavy quark energy loss and jet correlations in pQCD with the 
ones in AdS/CFT [HI gl, [H, M IH, El . 

In this work we study the propagation of a very energetic heavy quark through the SYM plasma at strong 't Hooft 



coupling and at temperature T. Our starting point is the trailing-string picture of Herzog et al 21| , which allows to 
compute the energy loss of the heavy quark in a static infinite-extend plasma. The heavy quark, in the fundamental 
representation, lives on a flavor brane close to the boundary with a string attached to it, hanging down to the horizon. 
The classical dynamics of the string is given by the Nambu-Goto action and from the string shape the energy loss 
can be determined. In [U a constant electric field is imposed on the flavor brane to drag the heavy quark at a 
constant speed v. The electric field feeds energy into the string and is adjusted to compensate the energy lost into the 
plasma. This calculation gives the leading order contribution to the energy loss when the mass of the heavy quark is 
extremely large. However, it is not such a realistic model. 

In this paper we study what happens if one turns off the electric field, meaning if one stops pulling on the heavy 
quark. Then the energy loss of the quark is not compensated anymore, and the string actually losses energy, the heavy 
quark decelerates. Since we are working in the limit of large mass M, the deceleration is slow and v ^ ^/XT'^/M can 
be treated as a perturbation. The deceleration modifies the shape of the string which in turns modifies the rate of 
energy loss by a correction of order \f\T jM. This adds a higher order correction to v and so on. This is the problem 

that we solve. The calculation is carried out up to second order in We explain the physical meaning of each 

correction which adds insights to the process of the heavy quark slowing down in a strongly coupled plasma. Finally, 
in this picture the trailing string slowly gets straighter to eventually become that of a static quark. Of course we do 
not except our picture to apply that far. What will happen is that the heavy quark will eventually thermalize and 
when that happens, the discussion should be modified. However our calculation allows to estimate the thermalization 
time. 

The plan of the paper in as follows. In Section II, we recall the string equations of motion derived from the 
Nambu-Goto action. In Section III, we discuss the possible ways the string gains and loses energy, and define the 
medium-induced energy loss which we are interested in. In Section IV, we compute the first and second order 
corrections to the trailing string picture, when one stops pulling on the heavy quark and it slows down. Finally 
Section V is devoted to an estimate of the thermalization time. Section VI concludes. 

II. THE EQUATIONS OF MOTION OF THE STRING 

We consider the large iVc, small gauge coupling gyu limit 

A^c ^ oo , gyA/ ^ , A = 5yjvf^c finite, (1) 

where the 't Hooft coupling A controls the TV = 4 SYM theory. The strong couping regime means A ^ 1, and then 
the equivalent string theory in AdS^ x space is weakly coupled and weakly curved: 

QYM < 1 ^ < 1 and A > 1 ^ i? > , (2) 

where gs is the string coupling, 1^ is the string length and R is the curvature radius of the AdS's and S'5 spaces. In 
this limit of small string coupling and large curvature radius, classical gravity is a good approximation of the string 
theory and we can neglect excited modes along the S'5 space. The metric corresponding to the SYM theory at finite 
temperature is the AdS black brane solution in 5 dimensions which can be written as 



^^■5^ ^ \ TTTTTT ~ u^f{u)dt'^ + {dx^ + dy^ + dz^) | , /(u) = 1 - ^ , (3) 



du2 



where uy^ — ixT is a black-hole horizon. The corresponding Hawking temperature T is the temperature of the SYM 
theory. The coordinate in the fifth dimension u — r/R^ has the dimension of momentum and the SYM theory lives 
formally on the boundary at it = 00. 
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In order to have one flavor of quarks in the fundamental representation in the field theory side, one introduces a 
D7 brane in the string theory side, whose embedding covers the part of the ^^5*5 space with u > Um = 27rm/-\/A, 
m being the lagrangian mass of the heavy quark^. The heavy quark and the color field it generates are dual to a 
string attached to the D7 brane at u = Um and hanging down to the horizon. Points on the string can be identified 
to quantum fluctuations in the heavy quark wave function with virtuality ~ u. Indeed, the quantum dynamics in the 
SYM theory is mapped onto classical dynamics in the 5th dimension. More precisely, the string dynamics is given by 
the Nambu-Goto action 



S = —To J drdaC = —To J drday/ — det gab 



(4) 



where (r, a) are the string world-sheet coordinates and Tq is the string tension. We define (r, a) as a map 
from the string world-sheet to the five dimensional space time, and introduce the following notations for derivatives: 
= drXf^ and X'f^ = daX'^. The determinant of the induced metric — det 9^;, = —g is given by 



- dct gab - {xf^X'^y - [X^X,,) (X'^X;) . 
It is useful to define the string canonical momentum densities as follows 

(^x-xi)x'^-ix'-xi)x^ 



dC 



'9 



(x^xi^ Xf, - (x^x^) 



dX' 



-9 



(5) 



(6) 



(7) 



Now let us consider an arbitrary variation 5X^^ in the action. After neglecting the surface terms, we can obtain the 
equations of motion of the string according to the variation principle 



5S 
'SX^ 



= 



dC 



d dC 



d dC 



dr dXt" da dX't" 



= 0. 



(8) 



When one chooses a static gauge by setting (r, a) = (i, u), and defines X^ = {t, u, x (t, u) , 0, 0) (we assume that the 
quark moves along the x direction), it is straightforward to find that 



-det gab = R'' [I 
The string canonical momentum densities reduce to 

— XX J 

-l-{x'fu^f 





u^f{u)x'^ 



-x' f 

-l + (i)2/-l 
X x' f 



(9) 



(10) 



For X^ = X {t, u) and — t, Eq. ([5]) gives two equations which represent the conservation of momentum and energy, 
respectively: 



dr da 
dr da 



0, 
0. 



For X^ = u, one obtains 



dr 



da 



-To- 



X- ul 



3 ^'2 



+ u-'x' 



(11) 

(12) 
(13) 



^ In fact, the relation between m and Um is modified at finite temperature. See Eq. (3.3) of Ref. [2lj. However, this modification will bring 
only higher order terms in our large mass expansion compared to the terms we calculate. 
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Due to the curvature in the fifth dimension, ^ generates a tidal force on the string in the u direction. However, 
among those three equations of motion above, there is only one independent equation due to the gauge fixing. One 
can just solve Eq. (fTT|) . and check that the solution x {t, u) is also the solution of the other two. Therefore, derived 
from Eq. (jlip . the equation of motion of the classical string in the x direction can be written as 

d fu^fiu)x'\ 1 5 , ^ >^o. (14) 



du \ \p-g ) f (u) dt 



III. THE ENERGY OF THE STRING DURING THE HEAVY QUARK DECELERATION 

The energy flow along the string at the position u = <j towards smaller u is given by 

u) = < = ^^ia;V/(u). (15) 
The energy density of the string is given by tt^ as: 

e{t, u) = -< = [1 + x'^u^fiu)] . (16) 



We consider a string which moves at a velocity v in AdS^,. Generally speaking, at finite temperature and velocity, 
there is always an event horizon in the fifth dimension which divides the trailing string into two parts. The upper part 
genuinely belongs to the dual of the heavy quark and its wave function, while the lower part is part of the plasma. 
This horizon is located at^ u — Ug where 

Us = ^/^Uh with 7 = ^ ■ (17) 
V 1 — V 

On the gauge theory side, this means that the part of the string above its corresponds to highly virtual fluctuations 
part of the heavy quark wave function while the part of the string below Us corresponds to longer-lived fluctuations 
which became emitted radiation. In the following we study the system formed by the upper part of the string, dual 
to the heavy quark dressed by the highly virtual fluctuations. Therefore, we consider the energy of upper part of the 
string 

E{t) = / du e(t, u) = - / du Til = Toi?4 / du i_' — . (18) 



Then, from (|T8l) . one calculates the energy change of the system during dt 

dE = d(^- J du TTl{t,u) 



Trl{t,Us) d{us) — dt du 9t7r^(i, u) 



— T:l{t,Us) d{us) + dt j duduT^ti^T^) 

J Us 

= TT]:{t,Us) d{Us) - TT^ {t,Us)dt + TT^ {t,Urn)dt , (19) 

where we have used Eq. p2)) to go from the second to the third line. Hence, the global energy conservation equation 
of the system during dt reads 

dE^5E\^^ -<^{t,Us) dt + ^t,u„^) dt, (20) 



^ This has been shown expUcitly in Ref. '43l. is also the location where the local speed of light in AdS space coincides with the 
propagation speed of the string. From the heavy quark point of view, the event horizon is at Us ■ 
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where 

S E\,^^^ ^ nl{t,u,) ^^^dt = e{t,u,) ^^^i^^dt (21) 

is the energy gain of the system during dt due to the increase of the size of our system. Indeed, Us decrease as v 
decrease. 

dt = dt (22) 

is the energy loss of the system during dt due to the energy flow along the string at m = it^, and 

^{t,u,^) dt = ^:"^{t,u,^) dt (23) 

is the energy gain of the system during dt due to the incoming energy flow at the top of the string at u = Um- 

In the heavy quark energy loss problem, these contributions can be interpreted as follows: $(i,u„i) dt is the work 
of external forces acting on the heavy quark and $(i,Us) dt is the energy radiated by the dressed heavy quark: 
fluctuations softer than the scale given by Us are freed into the plasma. As the heavy quark decelerates, Us becomes 
softer and more fluctuations are kept in the wave function. This leads to the energy gain 5E\^^ . 

In Ref. [2l| . an external force is acting on the heavy quark and adjusted in order to have a stationary solution, with 
V = Q and dE/dt — for the upper part of the string. Thus, the energy flow is uniform along the string, meaning 
$(t. Us) — '^{t, Um)- The same amount of energy is put into the string at the top than flows at the bottom into the 
horizon, or equivalently is radiatively lost by the heavy quark. However, this model does not describe the heavy quark 
energy loss so satisfactorily, since in reality the heavy quark will slow down when propagating in the plasma. 

In this paper we would like to consider the case where the external force is turned off and the heavy quark is allowed 
to decelerate slowly. In this case, we always keep <^{t^Um) vanishing, and we compute the energy loss. Although the 
total energy loss —dE/dt of the upper part of the string includes the contribution from the term S E\j^^ , we focus 
on the radiative energy loss $(t, m^) which is physically more important. For instance, this contribution could be 
compared to calculations performed in the pQCD framework [13, [HI [12, [3 • 



IV. PERTURBATIONS OF THE TRAILING STRING DUE TO DECELERATION 

A. The stationary solution as the leading contribution in a large quark mass expansion. 

Let us review briefly the main result of Ref. [2l|, , which is the stationary solution corresponding to a the heavy 
quark {um — > oo) pulled by a external force imposing a constant velocity. The stationarity condition gives the Ansatz 

X = xq + vt + vF{u) with i) = . (24) 

Since there is no time dependence in both x and x' , the time derivative gives vanishing result in the equation of 
motion. Thus, one can equate both the space-dependent and time-dependent parts of Eq. to zero, respectively. 
Therefore, one obtains 

'"^K/(Vi-|-f' ^''^ 

where Co is an integration constant. In order to maintain flnite and positive value of x'^, we should require that 
the numerator and denominator change sign at the same u value {u ^ Ug = ^/juh). Thus, this fixes the constant 

2 

Co = uf^vj, and yields x' = zfc ; where the -I- sign corresponds to a flow into the horizon at u^ and the — sign 
corresponds to a flow out of the horizon. By requiring the incoming boundary condition, we pick up the solution with 
the -|- sign and discard the other solution. 
Then, the shape of the string is found to be 

2uh 



TT 

2 



tan 



coth 



with vF'{u) 



vut 



(26) 
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and with the boundary condition F (u)\^^^, 
string and the energy flow along the string: 



0. Using that solution, one can compute the energy density of the 



e{u) = 



[1 + u'^f{u)x''- 



A 4 ' 



xx'u^ f{u) = TqR^ l^^u^ 



(27) 
(28) 



The energy flow along the string is constant, and the energy stored in the upper part of the string (u > Ug) is 

f du e{u) = ToR" 7 {"m - Us + v^ul[F{u„,) ~ F(u,)]} . (29) 

J Us 

The energy conservation equation (|^D|) is trivially satisfied: dE/dt — S /dt — and $(«,„.) = ^(ug)- 

When the external force is turned off, $(i,Mm) — and in order to keep ^20p satisfied, i) must be non zero (unless 
the mass of the heavy quark is stricly infinite). Indeed dE/dt is no longer zero: 



dE 
Itt 



7+1 



7^ + 1 
272 ' 



(30) 



and this imposes that in the large u,„ limit, —v oc u^/u^ This follows from the fact that in the energy conservation 
equation (PO)) . the finite term (when u„ 00) ^{t,Ug) — TqR^ jv'^uf^ + 0{v) in the right-hand side must be 
compensated by a finite piece in dE/dt. There is indeed no finite terms coming from the other contribution 



5E\ 



dus 



dt 



-e{t, Us 



— - — = —J^oR V V Us + 0(v ) 



(31) 



Therefore we now treat v as a, small parameter in the heavy quark deceleration problem. On the field theory side, 
this corresponds to a small (r\/A)/(2m) expansion according to the AdS/CFT dictionary. 



B. First order in solving the equation of motion 



Since the deceleration is very slow for heavy quarks, one can view the string as a quasi-static trailing string at 
every moment. Its shape should not be too different from the trailing string solution found above, except close to 
u — u„i where corrections are of order one, as is the case for $(<:, Um) which is now vanishing. Thus, we perform a 
perturbative expansion of the shape x(t,u) of the string as a series in the small parameter v/uh- In our quasi-static 
picture, we can assume that the coefficients of the expansion depend on time only through v{t). Hence, we write the 
generalization of the Ansatz including the first order correction as 



x{t,u)=xo+ / v{t')dt' + v{t)F{u) +v{t) C(u,v{t)) 



(32) 



After neglecting the and v terms, which are suppressed by one more factor of we reach the equation of motion 
at first order in — (the terms which are proportional to w). 



Solving Eq. dSS]) yields 



u^ulF{u)-v^ul^^C, 



1' [u-C] 



7'^u 

7R 



-I- F'{u) F{u) 



(33) 



(34) 



3 The 0{v'^) terms come from the fact that after switching off the external field, the energy of the string differs from II29I I by a contribution 
of order v. 
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As for the previous order, we determine the integration constant Ci by requiring a regular shape of the string at 
u = Ug, and thus Ci — Us- After the final integration and requiring that C(^)L=oo ~ ^' "-"^^ S^^® 



C{u,v) = 



2u2 



tan 



log 



U + Us 



F2 (u) 



(35) 



Strictly speaking the conditions C(^^)ljj=oo = ^ ^"^^ ^(^)L=oo = ^ imply that is not really the velocity of the 
quark, which is supposed to live at m = Um- However, as we show a 'posteriori in appendix, the difference is relevant 
only at higher orders in our large mass expansion. 

For any arbitrary point on the string, one can compute the energy flow <i>(t, u) which runs through it from the top 
part of the string to the bottom. First at u = Um, all the contributions in x' are of the same order, and one obtains 



■ ^ T.T. U J(U) 

ToR^ 



V {vul+ vCiUm:V)u'^f (Um)} 



Imposing the Neumann boundary condition <i>(t, Um) = on the flavor brane gives 



= ToR^-fv^ul 



1 



O 



The deceleration of the string at first order is then given by 



(36) 
(37) 

(38) 



J^{u.,ri ~ Us) 



(39) 



This equation allows us to specify our expansion parameter, it is actually Uh/{um — Us), and not Uh/um- In particular, 
our perturbative expansion breaks down if Us is too close to Um , meaning if the point which determines the energy 
loss is too close to the flavor brane'*. In order to interpret this parameter on the gauge theory side, let us come back 
to the energy of the upper part of the string. The first order contribution is 



E{t) = roi?2 ^(^^ 



i + o 



Uh 



(40) 



On the gauge theory side, this is the energy of the color field of the quark. It is smaller than in the vacuum due to 
the screening of soft fluctuations by the plasma. As in Ref. [2l|, , we shall use an abuse of language and denote 
E the thermal mass of the quark M. It is different from the Lagrangian mass m because of the flnite temperature, 
but the interpretation as a quark mass on the field theory side is not entirely clear. In particular, while for the heavy 
quark at rest all three definitions of the mass M coincide, for a moving quark the 7 dependence proposed in (21I . l49j 
is different from ours. In our case, the 7 dependence accounts for the fact that as the heavy quark decelerates, softer 
modes are allowed (not screened by the plasma anymore) in his wavcfunction: 



Us) ■ 



(41) 



It is also defined so that —Md-f/dt — ^{t,Us), the radiative energy loss, which for now can be checked at leading 
order. In terms of M, we can now write our expansion parameter as {T^/X) / (2M). In appendix, we give the second 
order and third order (this one is obtained from the correction C(u, v) to the trailing string) corrections to the energy 
and therefore to the mass. 

Now we can compute the energy fiow for u <C Um, using equations ((32|) and (|34p : 



$(t, u) 



ToR' 



Xx'u f{u) — TqR "fV Uf^ 



V 



l + (l+7^)^F(.) + '^'^";"-^ 0. , 
V vui \ui 



(42) 



* For u < Um, this string picture is modified anyway [33. l47l. |43| . The string is actually a narrow tube: the interaction of the string with 
the D7 brane deforms the brane in the neighborhood of the endpoint of the string over a length scale of order 1/M. 
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By evaluating equation (|42p at u = Ug and inserting equation (|39p , one finds the outgoing energy flow from the upper 
part of the string 



<^>{t,Us)=ToR^-fv^ul 



1 



O 



Using formulae (pi]) and it is easy to see that the energy conservation equation 

dE _ . SE\ 



I dus 



dt 



(43) 



(44) 



is verified, and to recover the net energy loss of the system at first order ((30l) . Using the AdS/CFT dictionnary, one 
deduces the energy radiated as soft gluons by a dressed heavy quark moving through the SYM plasma at a speed v 



2,71 

~ (ttT) 

7>1 




(45) 
(46) 



As F{u) is negative, the first order correction to the soft gluon emission $(t,Us) is always positive, meaning that 
the rate of energy loss is stronger for lighter quarks. Interestingly enough, in the ultrarelativistic case 7 3> 1, the first 
order correction to $(t,Ms) is further suppressed by a factor 7""^/^, which is unexpected. In the non-relativistic case 
7 > 1, due to the expansion F{u) ~ log((u — Uh)/uh) when u Uh, the first correction to ^{t, Us) is dominated 
by the logarithmically enhanced term. It shows that our large mass perturbative expansion breaks down at very small 
velocities. And even if due to the prefactor that logarithmic singularity does not lead to a divergent energy loss 
in the small velocity limit, but to a vanishing one as expected, in reality 7 never reaches 1. The heavy quark will 
thermalize before this happens and already then the picture has to be modified. 



C. Second order in solving the equation of motion 

At second order in Uh/{um — "s), we assume the solution of the string to be of the following form 

x{t,u) =xa+ [ v{t')dt' + v(t)F{u) + v{t)((u,v(t)) + v'^{t)xi {u,v{t)) + i{t)x2 (u,v{t)) 



Thus one finds 



= v'^ + 2vvF{u) + 2vvC{u,v) + v'^F^{u) + 2vv 



dv 



a;'2 = v^F'^iu) + 2i)vF'{u)C{u, v) + 2vij^ F' {u)x'i {u, v) + 2vvF' {u)x'2{u, v) + v^C{u, v) + 
Expanding the time derivative part of the equation of motion up to second order, one finds 



1 d 



f{u) dt 



1 



/(") 
■2 



7 f + jv 



v'y^F{u) 



fin) 



/(") fin) 



fin) 
+ 1 - 



■ 2^n^ 

fin) 



{2 + 3-f^v^)Ciu,v) + v 



dC{u,v) 



dv 



(47) 



(48) 
(49) 



(50) 
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Expanding the spatial derivative part yields 
d f u'^f (u) x' 



du 



= ivdu 



iti \ 2 2 4 a// \ , 2 2 2 ^(^) 



7H 



( Fju) 



7 vu,^ 



IF^u) 1 



2 f{u) 2 



f (u) dv 



fin) 



\C{n,v) 



(51) 



Comparing the coefficients of and v gives the equations of motion for xi a-nd X2, respectively^: 
u^f{u)x[{u, v) - j'v'utx'iin, v) + j\u^f{u)C{u, v) ( ~ ulC{u, v) 



+du 



~1W 



lF^{u) 1 

— —I 

2 f{u) 2 
2 

fin) ^ f{u) 



V dC(u,v 



3 ,1 q 9 I F{u) 2 /■!( ^ 



9 9 



dC'{u,v) 



dv 



and 



4 r/ \ I / \ 224// \i 222 Ci^) 



fiu) 



F(«)+7V^^"^ -2.,2.,2^// 



fiu) 



7 V u^^Ciu,v) 



(52) 



(53) 



First, let us integrate over u once on both sides of the equations. For Xi('u, v) and X2(''^i '^^^ g^ts, respectively: 
{u^ - ut) xU«, v) = -j'vu^f{u)C'{u, v) - n^X'iu, v) 

\2 f{u) 2 f{u) dv J 



"^ 4 3 2 ( F{u) 2*// \ \ 2 2 Ff) I o 2 21 

- 2'T « V ^ - f7 [2 + 37 f ' 



C(^^,z;)--^^2 



-"f'^vu\ [2 + 37^v^] (m, w) + 3w7"'t7 {u) — v'^^'^cj) {u, v) + C2 , 



and 



/'4 4\// \ 222 '^Z I 2 / \ 

-u,)X2[u,v) = -J V U^—— + a {u) 

-v^uW \(:m~]^f^ {u) 



- I'^v'^ulg (u, v) + C3 



where where C2 and C3 are integration constants. We have also introduced the following functions 

1 u — 



g {u, v) — Uf^ du 



du 



F ju) _ ul 
'fin) 2 



F^ (u) + uF (m) + - In 



1 . u'^ + ul 



(j) {u, v) = u\ I du 



1 dC'M 

'f{u) dv 



l^h I du 



4 — u\ 
1 Au^ + iu^Us + 2ulu + ul 
l{u) 2{u + u,f {u-^ + ulf ' 



(54) 

(55) 
(56) 

(57) 
(58) 
(59) 



^ Strictly speaking, there is physically only one unknown function in our Ansatz l|47| l at second order, vxi + ^Xi, but we choose to 
determine xi and X2 separately. 
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with 



giu,v) 



Uh 



Uh 



2ut tan 



2u)j tan 



+ uj. (uh - Us) log (U - Uh) 



4u2 - uj) 



[ul (uh + Us) log {u + Uh)~ 2ul log (u + Us)] 



Uh 



4«2 (^,4 _ ^4) 



J- [~uluh log {U^ + ul) + ul log + u2)] 



(60) 
(61) 
(62) 



By requiring that Xi{u) and X2(") are finite at u = u^, one can determine the integration constants. They are such 
that in the large u limit for the right-hand side of ([54|) and ((56|) . one finds 



{u''-ut)x[{u,v) = 5- 



37^ 



O 



{u^-ut)Mu.v) = C3+(l-^)4?^+(^(^) + 



with 



1 57r 5 log 2 
4 ^ "8 4~ 

TT log 2 

4 2~ 



■i;7 



7 



(63) 
(64) 

(65) 
(66) 



We shall not integrate the equations of motion any further to get the shape of the string, the knowledge of x'i{u) 
and X2(w) is enough to compute the energy flow $(t,M). First at w = u„i, all the contributions in x' are of the same 
order, and one obtains 



$(i,Mm) 



J-9 

Vi — t'' 



xx'u f{u] 



■.V {vul + [v('{u„i, v) + V^x'l ("m, v) + {Um, v)] U^f [u^)] 



Imposing a vanishing flow at the top of the string gives the second order correction to i) : 



vul 



1 TT log 2 
72(u„-Us) ' 7(Mm-Ms)^\4 



(67) 
(68) 

(69) 



Note that in the second order correction, we only gave the large 7 result®, in order to get a lighter expression. We 
will also need the second order expressions 



„-.2 



v'il-v')ut 



v{l-3v')ut 



(70) 



7^(Wm-"s)^ ' -f^iUm-Us)^ 

Now we can compute the energy flow for u ^ Um up to second order in Uh/ {um ~ Us), by expanding y^—g, x, and 



While formula l|17| l is valid for the trailing string solution, the first order correction to the string shape modifies the 7 dependence of 
Ms , the point where the local speed of light in AdS space coincides with the propagation speed of the string. In principle this should be 
accounted for in the second order calculation, however in the large 7 limit, these corrections are always subdominant. 
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x'u fiu) up to second order: 



TqR'^ . I 4 



-9 



xx'u f{u) 



(71) 



= ToR^ttT) 



2 W 



+ToR'^ (ttT) 



2 V 



9 4 4 

7 u — u\ 



2 



9 4 4 
7 W — 



V 

v(^{u, v) 

V 



vut 



2 V 



2 



• 2 2 

u 7 



C'^(u,w)u4/(u) 



2 W 



2 



2/ (m) 2 



2 W 



1)^7^ [F{u) + 

^ 9 2 



/(«) 



vul V dv — ui 



(72) 



The energy flow at u = Us does not depend on x[ (u) or X2 ^'^^ found to be (in the large 7 hmit where the 
second order correction to v is subdominant) 



$(i,u,)^Toi?2 uIjv^ 1 + 
where we have used the results 



Uh 



87^/2 Urn - Us 



dv 



log 2 7 

8 4 ^ 32 ' ' 



o 



1 /I TT 1 



'"^ 2ul\4 ■ 4 2 



log 2 



(73) 



(74) 



In order to rewrite this in terms of the heavy quark mass, and other gauge theory quantities, one should first 
determine M up to second order. From formula (j29[) . one finds 



M = ^{Um - Us) (1 

Zn 



Um Us 



F{Us) ] J^^-{^{Um - Us) ( 1 



73>1 



2tt 



v'^Uh 



il^^'^iUm - Us) 



(75) 



which is again consistent with —Md'j/dt = ^{t,Us). Then in formula ([73|) . when introducing {um — "s) ^ = 
VX/{2ttM) [l + 0{uh/{um-u s))] , one sees that what was first order in Uh/{um — "s) will introduce a second order 
correction in terms of ^fXT /{2M). However, in the large 7 limit, this contribution is subdominant. Therefore one gets 
for the radiative energy loss 



<^{t,Us)^^^{T:TfW 

7»1 



1 Vat (tt 1, ^ 7\ (\T^\ ,^/A='/2r3 



373/2 2M 



- ^4^°^'- 32] niM^j 



O 



SAP 



(76) 



The second order correction is also positive (f — -j log 2 — ~ 0.00066) but very tiny. This continues to enhance 
the energy loss for heavy quark with lighter mass. Also, one can check that the second order correction is always 
small as long as Ug = ^/jUh <^ Um- It is also interesting to notice that the second order correction is proportional to 

-^a^ [sol-lsH when a is identified as "f'^v in relativity. One can interpret this term energy lost by radiation due to the 
deceleration. 



V. THERMALIZATION TIME 



To finish, let us estimate the thermalization time of the heavy quark, and show that our perturbative expansion 
can be trusted, provided is large enough and the heavy quark slows down in a quasistatic way. In order to show 
how our calculation improves the previous estimations, let us rewrite formula (|69p in the following way: 



dv 



-v{l~v^) 



1 



7 



1 



log 2 



(77) 



12 



fi = 10 and = 3 




1 r 

1st order 
2nd order 
eq. (78) 



60 



70 



FIG. 1: Numerical solutions of equation (|77|) without the square bracket term (first order), with the square bracket term 
(second order), and the solution (|78|l proposed in [2ll. [2^. Our results indicate a faster heavy quark thermalization. 



where i ~ ttT t and /i = u„i/uh = 2m/ {^/\T) is expressed in terms of the bare mass of the heavy quark. Note that 
/i ^ in order for our perturbative expansion to be vaUd. 



In Ref. [2l|, l22|, the equation proposed for v is fi dv/dt — —viX—v^) and the solution is 



v[t) 



\.§ + (l-.§)exp(^) 



(78) 



with = \/ 1 — 1 / 7g the initial velocity. This naturally implies the thermalization time scale to be . However 

the equation proposed in Ref. [Ill, [l^ is different from our first order equation {^—^J^) dv/di = — w(l — w^), which 
indicates that the thermalization time is rather given in terms of the thermal mass of the quark: 



2M 



(79) 



Strictly speaking M slightly varies with t, from m — T\/\^q/{2t:) to about m — T^/X/ {2tt), which bounds the ther- 
malization time. In any case it is a smaller time than ^"^^ , and this leads to a faster thermalization. 

Let us point out that this difference comes from the fact that we had to pick a particular point on the string to 
evaluate the energy loss, while in the stationary problem solved in Ref. [U, [l^], it doesn't matter at what value of u 
the energy flow is computed because the flow is constant when w = 0. In our study, we actually let the heavy quark 
slow down and the point on the string where the energy loss is evaluated becomes crucial. Our choice to evaluate the 
radiative energy loss at u — Ug is motivated by the fact that the part of the string below Us is not causally connected 
to the part of the string above [43], and therefore corresponds to emitted radiation on the gauge theory side. 

In Fig.l, equation (|77p is solved numerically at first order (without the square bracket term) and at second order. 
The comparison with the solution ([75)1 confirms that our calculations imply a faster thermalization: the time it takes 
to reduce the initial speed in half is about 2/i/ (ttT), instead of 3^/{ttT). We also see that the second order correction 
is quite small, although one should keep in mind that the dashed curve should be modified for small values of 7 : the 
second order correction in (|77p is only valid in the limit 7 3> 1. 



VI. CONCLUSION 



In this paper, we study the radiative energy loss of a slowly decelerating heavy quark moving through a strongly 
coupled SYM plasma by using AdS/CFT correspondence. The calculation is done in terms of perturbation in ^/XT/M, 
where A is the strong coupling of the SYM theory, T is the temperature of the plasma and M is the thermal mass 
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of the heavy quark. The shape of the string traihng behind the heavy quark in the fifth dimension determines the 
medium-induced energy loss of the quark. At zeroth order, the heavy quark propagates at a constant speed, and 
the traihng string l(24 |l - ((26|l is that found in Ref. [Ill, [l^. We determine the corrections to the shape of the string 
up to first order, formulae ([5^ . ([55)1 and ([5^ . and the radiative energy loss of the heavy quark up to second order 
(|73p . In this case we only gave the large 7 result explicitely, but the exact expression was also found. The physical 
interpretation of each correction is provided. Especially, the second order correction can be understood as energy loss 
due to the decelaration. Last but not least, these calculations allow us to evaluate the thermalization time of the 
heavy quark moving in the strongly-coupled SYM plasma. Our parametric estimate is given by d75|) . 
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APPENDIX A: VELOCITY OF THE QUARK vs. v{t) 

Using our Ansatz (|T7)) . the velocity of the heavy quark, which is the velocity of the string at the flavor brane, is 

i{t, u^) = v{t) + vit)F{u„,) + i{t)C{u^, v) + v^t)d.,C{u,n, v) + ... . (Al) 
On the other hand, we have the expansions 

FM = -^ + o(4-) , (A2) 

aum,v) = -J^ + o(^) . (A3) 

Therefore, the velocity of the quark reads 

x{t,Urr.)^v{t)+0(^^^ . (A4) 

The discrepancy between v(t) and the velocity of the quark appears only at higher orders in the large mass expansion 
than the ones considered in this study. 



APPENDIX B: CHECKING ENERGY CONSERVATION UP TO SECOND ORDER 



According to Eq. ([TC)) the energy of the system is 



E = 



du e{t, u) = TqR'^ 



du 



1 + u^f{u)x'^ 



Inserting Eq. (|47| . one finds 



= 1 



1 



4 + vvY 1 + 



F{u) 
fin) 



- uiCiu, 7) + 2vvuiCiu, v) + 0{ — 



and thus, neglecting terms suppressed by powers of Um, one gets 



E = ToR'-f u,„ -Us- v'uiF{us) - vi - (7' - 2)<C(«., t-) + r<j{us) + {Y - l)^F\us) 



(Bl) 



(B2) 



(B3) 
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Using 



SE\ 



dus 



dt 



Us^fVV 111 

-e(t,Us) — - — =ToR V 



(7^ 



1) 



27 u. 

1 . „4 



4m., 



O 



{Um - UsY 



(B4) 



and substracting (formula (|73p '). one can check that the energy conservation equation is verified, meaning 

that one recovers the net energy loss dE / dt from derivating (|B3|) . 
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